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ABSTRACT 

This document develops an analytical/numerical approach to identifying and 
correcting the aberrations introduced by a general displacement of the feed 
from the focal point of a single offset-paraboloid antenna used in deployable 
radiometer systems. A IS^neter reflector with 18*meter focal length is 
assumed for the analysis, which considers far-field radiation pattern quality, 
focal-region fields, and aberrations appearing in the aperture plane. The 
latter are obtained by ray-tradng in the transmit mode and are expressed in 
terms of optical notation. Attention is given to the physical restraints im- 
posed on corrective elements by real microwave systems and to the inter- 
mediate near-field aspects of the problem in three dimensions. The subject 
of wavefronts and caustics in the receive mode is introduced for comparative 
purposes. Several specific examples are given for aberration reduction at 
eight beamwidths of scan at a frequency of 1 .414 GHz. 
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unit nonnal to a surface 
incident ray vector (in context) 
reflected ray vector 
optical invariants 

transverse image displacement parameter, radial and azimuthal variables 
phase perturbation (in context of aberrations) 
phase coefficient (in context) 
aberration indices 

Cartesian coordinates of fimdamental reference frame 
beam squint angle 
beam deviation factor (in context) 
aperture diameter 

involute, evolute of differential geometry 
wavefront 

caustic (surface, line, or point in context) 

radius of Airy disc 

feed displacement in focal plane 

radial and azimuthal variables of offset paraboloid 

focus or focal length (in context) 

parameters of elliptical beam 

a displacement vector 

Euler angles 

polar and azimuthal angles for elliptical-beam feed 
decibel 

unit Cartesian basis vectors 
generic symbol for surface 
wavelength 

vector sheet cxirrent of physical-optics theory 
incident magnetic fleld 
polarization moments 

scalar weighting factor of elliptical beam (in context) 
phase argument 
qiiatemion rotation operator 
degenerate quaternion (vector Wj ) 
rotation angle in Euclidian space (Ej ) 
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THE CORRECTION OF ABERRATIONS COMPUTED IN THE 
APERTURE PLANE OF MULTIFREQUENCY 
MICROWAVE RADIOMETER ANTENNAS 


INTRODUCTION 

The relative success of an attempt to compensate for the undesirable effects introduced by the dis- 
placement of an antenna feed from the focal point of the reflector may be evaluated by calcvilating 
(or measiiring) the far-field radiation pattern. Because this pattern is the Fourier transform of an 
aperture distribution, an analysis of the amplitude and phase in the aperture plane for a specified 
polarization may reveal most of the features of the far-fleld pattern by inspection. An aperture 
plane analysis is economical since integration is not required, it provides a large amount of informa- 
tion, and it leads to an early diagnosis of the aberration problem. 

In the present doaunent, the transmit point of view is taken, and phase departures from an apertiuc 
plane of a paraboloidal reflector are determined and classified by analogy to an optical approach. 
The latter approach takes the receive point of view at a Gaussian focus and regards phase departures 
from a Gaussian sphere. The receive point of view is retained in the present dociunent to locate dis- 
placed focal regions and determine their size. This is accomplished by several means, including a 
simple three-ray trace method, a multiray caustic analysis, and a focal region mapping based on 
physical optics theory. It is shown that the caustic may, alternatively, be obtained as the evolute 
of the received wavefront. Subsequently, the placement and size of corrective feed elements is 
predicated on this knowledge of the focal region. 

ABERRATIONS 

The aberrations observed in the focal plane of an offset paraboloid in vacuo may be due to the wave- 
front associated with the feed, the position of that feed with respect to the focal point of the re- 
flector, or imperfections in the reflector itself (Ref. 1 , p. 1 86). Assuming that the feed wavefront 
is spherical and the reflector is an ideal paraboloid, the feed displacement alone is responsible for 
the aberration. In any event, a ray trace that obeys Snell’s law of reflection and accounts for optical 
path length yields a phase distribution in the aperture plane. 

The equations, 


and 


n X S, » n X S, 


n • Sj = -n * 


( 1 ) 


( 2 ) 


are fully descriptive regarding the angular relationship between the incident rays, reflected rays, and 
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surface normal, but the equation, 


(3) 


Sj = Sj - 2n(Sj ' ID 

is probably more useful. (See Ref. 2, pp. 354-355.) 

A convenient notation for classifying aberrations in paraxial optics is one that employs the optical 
invariants of a symmetrical optical system, 

Ij = h*, Ij = p’, I 3 = hpcos 0 (4) 

and the expansion for the optical path difference 


OPD = a = ^ cos" « . (5) 

Ion ,11 

Here, h relates to lateral image displacement and p and 0 are radial and azimuthal variables in the 
aperture. See Reference 3, pages 46 and 47, Reference 4, and Reference 5, page 859 for details 
of the classification. Large aberrations lead to all orders in the expansion. Apodization and aper- 
ture blockage have also been treated (Ref. 3, p. 95). 

If certain constant terms, a^ > b^h^ , c^h^, and d^h^ , are admitted as “aberrations,” the classical 
aberrations of Gauss, Seidel, and Schwarzschild may be generated by simple factors, powers of 
the optical invariants. The aberrations of degree zero, two, four, and six emerge as products of 
these invariants with mixed degree of zero, one, two, and three. (See Appendix A of this docu- 
ment.) 

An Infra-Red Handbook (Ref. 4) defines ^ ^ equal to zero. This is explained as follows in 
Reference 6 , pages 87 and 88 . “It is customary in Hamiltonian theory to postulate that in the 
wavefront expansion the reference sphere for each field angle has its centre at the Gaussian image 
point.” Also, “The constant term is omitted because in the definition of wavefront aberration we 
have assumed that both the wavefront and the reference sphere are chosen to pass through the 
centre of the pupil.” In the present aperture plane analysis, the Gauss “aberrations” are retained. 
This was also done in Reference 5, page 858, and Reference 7, page 18, where three, six, and ten 
aberration terms are associated with Gauss, Seidel, and Schwarzschild, respectively. 

The two views regarding the optical path difference (OPD) are illustrated in Figure 1 for an optical 
system and an offset paraboloid microwave reflector. It is noted that the offset reflector is in- 
herently an asymmetrical geometry, but is a part of a parent paraboloid that is a symmetrical con- 
figuration. The development reverts to the latter, on occasion, and then proceeds after invoking a 
truncation argument. 

The primary aberrations of Seidel are illustrated in Reference 8 , page 212, as surfaces lying above 
or below a drcidar disk representative of the aperture plane. (They are not reproduced here.) It is 
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noted that the Cartesian aperture coordinates are 


X = p cos 0, y = p sin 0 (6) 

so that aberrations such as the coma term may be rewritten to advantage as 

to exhibit cross-gradients among the phase variations. 

CAUSTIC SURFACES AND WAVEFRONTS 

The receive point of view may be taken to obtain an appreciation of the behavior of the offset 
paraboloid for various angles of wave arrival, initially in the limit of zero wavelength, by means of 
caustic surface analysis. Detailed discussions of point, line, and surface caustics are found in Ref- 
erence 8, page 792; Reference 9, page 498; and Reference 10, page S26. “One of the earliest dis- 
coveries in optics (F. Maurolycus, 1S7S) was that the rays of a normal system are tangential to a 
surface, the so-called caustic surface, and that the general ray is tangential to the caustic in two 
points.” (See Ref. 9, p. 156.) The caustic has physical reality (Ref. 1 1 , p. 79). It may be defined 
as the envelope to a system of rays of an orthotomic system or the locxis of the principal centers 
of curvatrire of a wavefront (Ref. 12). Used in conjunction with a “spot diagram,” the intersection 
of rays with a chosen image plane, the caustic sxuface approach leads to useful impressions of geo- 
metrical image quality and intensity (Ref. 6, p. 93). 

The offset paraboloid geometry of the present discussion is shown as Figure 2, together with the 
surface-generating equations and assumed parameters. A system of three rays at plane-wave inci- 
dence 0** and >8** is shown in Figure 3. From this figure, it may be inferred that the focal axis and 
focal plane of an offset paraboloid are not the same as for the parent paraboloid, although the 
geometrical focal point is the same. Furthermore, the central ray of the three-ray approach is seen 
to be strongly indicative of the approximate focal axis (AFA) for various angles of wave arrival. 

The approximate focal plane (AFP) was, by analogy to parent paraboloids, initially taken to be 
orthogonal to AFA and passed through the focal point (F). The AFP was not obvious for general 
angles of arrival since the displaced pseudofocus (F') is not known prior to a focal region mapping. 

A ray trace based on 1 S rays (Figure 4) maps an arc of the caustic surface in the x-z plane for a 
plane-wave an^e of arrival of (-8*^). An approximation formula relating the angle of arrival (d,) to 
the lateral feed displacement (A) of a parent system of focal length (F) is: 

KA 

sin fl, » — . (8) 

Here, K is a beam-deviation factor (Ref. 1 , p. 488). Further, 

0.92 < K < 1.00 (9) 
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for 


0.50 < F/D <1.00 (10) 

of the parent system. If the approximation formula is applied to the offset geometry, and A is now 
taken to lie in the inclined focal plane AFP(0“) as shown, the location of the pseudofocus (F') may 
be estimated. It is seen to fall on the caustic arc for = -8°, but is not in the center of that arc. 
Rather, it is situated near one end of the caustic arc (K « 1 .0 assumed, F/Dp = 0.5, F/D^ = 1 .2). 

Retximing to the parent paraboloid, it is instructive to further develop the caustic, here based on 1 6 
rays, and identify the bounding rays (1 , 7, and 16) for the offset and parent surfaces as seen in the 
x-z plane of Figure 5. The caustic arc is nearly circular for the truncation of the parent paraboloid 
shown. This fact may be exploited to develop the involute (INV) of the caustic (w), or reflected 
wavefront (A), since parameterizations by arc length are sometimes difflciilt to formulate. The 
caustic (cj) is the evolute (EV) of wavefront (A). These concepts will be developed subsequently 
for the offset paraboloid geometry. 

For the present it is merely noted that: 


EV(A) = w 


( 11 ) 


and 


INV(w) = A 


( 12 ) 


It follows that: 


INV (EV(A)J = A 


(13) 


and 


EV [INV(cj)] = cj (14) 

where u may have branches (cu^ and cj,) as mentioned previously. 

The nearly circular caustic arc of Figure S is actually part of an inclined and translated curve known 
as Tschimhausen’s cubic (Ref. 2, p. 25). It has been associated with corrective lens methods for off- 
set paraboloids. Figure 6 shows the analytical caustic curve, which is equivalent to the envelope of 
rays tangent to that curve. The offset surface and incident plane wave comprise an asymmetrical 
geometry, but Tschimhausen’s cubic is a symmetrical function given by: 

(x^= y"(y"-b)2/3b (15) 

b = 9a sin . (16) 
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The inclination of the (x", y") coordinate frame equals d, and the translation follows from 


(y')^(y'-9asinfl) 
+ acos0 

(27a sin 8 )'^ 


( 17 ) 


when y' and y" are indistinguishable, upon setting (y') equal to zero. (See References 13 and 14.) 
The rotation transformation is given by: 


X = x'cos® + y' sihfl (18) 

y a y'cosd • x'sihd . (19) 


Finally, to arrive at the coordinates of Figure 6, the correspondence. 


xof(18)-a-»z (20) 

yof(19)-*-x (21) 

is made. It can be seen that the appropriate segment of arc on the Tschimhausen cubic curve is 
determined by the parent or offset paraboloid truncation. For the parent and offset reflectors of 
Figure 6, the bounding tangents intercept the cubic curve in the essentially circular part of that 
curve as shown previously by the envelope of the caustic (Figure 5). The cubic is extended here to 
display its symmetric form, principal axis intercepts, etc. 


Returning to the offset paraboloid and the microwave aberration problem, a set of caustic arcs is 
now developed corresponding to plane-wave angles of incidence equal to 0", -4®, -8®, -16®, and 
-32®. Figure 7 shows the bounding and central tangent rays for each angle of incidence and the 
associated arc along each cubic curve. The shaded zone illustrates the continuous and rapid increase 
in the caustic arc length with increasing (absolute) angle of incidence. An aberration problem at- 
tending the larger angles of wave arrival may therefore be represented by the phase departure from 
an aperture plane (ray-trace program data), the equivalent expansion (equation 5), or a part of a 
caustic surface. For simplicity, caustic arcs were illustrated here. A detailed development of the 
caustic surface for the offret paraboloid will be deferred to another Technical Memorandum. An 
appreciation of the complexity of the undertaking, which includes phase information on the caustic 
surface, may be obtained from Reference IS. 


FOCAL REGION MAPPING 

In an attempt to determine the focus of the offset paraboloid for -8.7® plane-wave incidence sev- 
eral methods, including complex-vector physical optics theory, were employed. The focal-region 
electric field was first computed for axial plane-wave incidence to verify the inclination angles of 
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the focal axis and focal plane. Figure 8 shows agreement with previous ray-trace results. A com- 
parison of the radius of the Airy disk for the offset paraboloid, surface integration sampling interval 
equal to 2.0 X, is made with the approximation formula for parent paraboloids: 



1.22 FX 
D 


(122) (18.0) X 
(15.0) 


1 .46 X 0.31 m (theoretical) 


( 22 ) 


« 1.75 X« 037 m (from Figure 8) . (23) 


The physical optics formulation and the parameters used to compute the focal-region fields at 0.0° 
and -8.7° plane-wave incidence are given in Appendix B of this document. 

Figure 9 shows the result of the -8.7° plane-wave computation, surface integration sampling interval 
equal to 1 .0 X, at 1 .414 GHz (X = 0.21 m). The maximum value of the electric field appears to be 
at (Xj , y^ , z^) = (2.8 m, 0.0 m, -1 .3 m), and the isophotes or constant-intensity contours are rather 
symmetrically disposed about that point out to 1 meter. Beyond 1 meter, the asymmetries of the 
offset reflector/inddent wave system can be seen. The radius of Airy disk for -8.7° plane-wave ind- 
dence has been greatly affected: 

R^(-8.7°) > 7Xor> 1.5 m . (24) 

Note that the caustic arc length at -8.0° inddence is 4.0 m (Figure 7). 

ESTIMATES OP THE DISPLACED FOCUS 

An estimate of the position of a displaced focus of an offset paraboloid for a given plane-wave ind- 
dence angle may be obtained by various means. A comparison of results is now made for the geom- 
etry of Figure 2 for an angle of plane-wave arrival equal to -8.7°. Figure 10 illustrates that the 
reflector and feed position information (x^, y^, z^) is obtained by several means. Figure 1 1 is an 
expanded-scale version of Figure 1 0. 

• y. . ^ = (2.56, 0.0, -1 .6) m 
A^Fsinfl/K;6 = 8.7°, F= 18.0 m, AFP(0°) 

F/D^ = 18.0/15.0 = 1.2, F/Dp = 18.0/36.0 = 0.5 
K ^ 1.0, 0.9, 0.8 (Ref. 1 , p. 488) 

Comment: valid for small d. 

(Fj) (X,, y,, z^) = (2.80, 0.0, -1.8) m, LI = 2.0 (X) 

Focal Region Map (Physical Optics) 

Comment: Generally reliable for larger d, but integration interval is large. Based on 
field symmetry. 
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(Fj ) (Xj > y* . ® (parabolic arc information only) 

o ^ 

z = ^-a\ > (0®) inclination 

* 4 (2F) 

Parabolic arc with vertex at origin (0). 

Attributed to P. Potter 

Comment: Arc has negligible departure from AFP (0®). 

4 

(F;> (x,,y,, 2 ,) = (2.8,0.0,-1.3)m,U = 2.0(X) 

(x^ , y, , z, ) = (2.8, 0.0, -1 .3) m, LI = 1 .0 (X) 

Focal Region Map (Physical Optics) 

Comment: Generally reliable for larger 6, two integration intervals used. Based on 
maximum electric field (E„.. ) value. 

max 

(F;> (x,,y,,z^) = (2.25, 0.0, -2.00) m 

Intersection of two tangent rays R,HRj (Figure 10) 

Comment: Experience has shown that the R^OR^ approach gives good results for 
small 8. 

(F') (x^,y.,z.) = (circular arc information only) 

R = F about vertex (V) as center (offset paraboloid) 

R - OC about midpoint (C) as center (parent paraboloid) 

Attributed to Silver, adaptation by author. 

Comment: Valid for small 8. 

(F;) (x,,y,,z,) = (1.8,0.0,-2.5)m 

Based on center of caustic arc bounded by R, , Rj 
(Figure 7, -8.7® interpolated) 

Comment: Large departure from other methods. 

(Fg ) (x, , y, , z^ ) = (not done) 

Focal Region Map 

Comment: Generally reliable for larger 8. Based on flatness of field (wavefront). 

(F' ) (x,, y, , z^ ) = (not done) 

Focal Region Map 

Comment: Generally very reliable for larger 8. Based on time-average Poynting vector. 

ABERRATION CORRECTION (I) 

On the premise that a given reflector antenna is able to concentrate the energy of an off-axis inci- 
dent plane wave in a rational “focal region,” it is necessary to move the primary feed to that region 
(Ref. 5, p. 860). Subsequently, it may be possible to adapt to the received field distribution and in 
some sense attain a “matching” feed that favors antenna gain, beamwidth, sidelobe levels, gain-beam- 
width product, a polarization characteristic, etc. Situations arise in which the reflector should be 
replaced with an adaptive array. It is assumed here that the reflector still focuses satisfactorily, but 
that one or more aberrations should be reduced, if not altogether eliminated, by combination 
methods if necessary. 
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Aberrations that are commonly encountered include transverse focus (Wj j ^ ), axial focus 

Ij ), coma (Wj 3 , Ij I3 ), and spherical ). See Appendix A-2 for additional aberrations. 

Transverse focus (Wj ^ j I3 ) is essential to beam squinting, is induced by a nonzero value of (h) 
representing lateral feed displacement, and is actually not an aberration affecting image quality. By 
the shift theorem of Fourier transform theory, it is shown that this phase gradient displaces a beam 
in space, but does not distort the beam. (See Ref. 1 6, p. 282, and subsequent discussion in this 
document.) The symmetrical aberrations, axial and spherical, may be dealt with by feed displace- 
ment. Nonvanishing lateral displacement (h) also introduces the asymmetrical comatic aberration 
^131 ^2 ^3 manifests itself in radiation patterns as raised and lowered sidelobes in the cut 

of the cubic phase gradient (Ref. S, p. 8S9). It constitutes a significant limitation for radiometer 
antennas. 

In an initial attempt to minimize coma, the feed function 


cos^9 

sin* 0COS* ^ + sin* 9sin*^ +cos* 0 

(1-e*) J 


(25) 


was developed and used to mathematically simulate five differentially phased, linearly polarized, 
elliptical beams with common phase centers and 3-dB crossover levels. The angular beam separa- 
tion was set at 9 degrees (A0) at 1 .414 GHz, spatial rotation being obtained by means of Euler 
angles (a = 90“, 0, 7 = 0®). The third primary beam (i - 3) was directed at the center (C) of the off- 
set paraboloidal reflector (Figures 12 and 13) as the pseudofocus (F') moved away from the focus 
(F) with increased plane-wave angle of arrival (main-beam scan angle). 


Since the electric field at theta (6) equal to 0® is unity, using equation 25, 


dB » 20Nlog 
10 


1 r 

■ + 10 log sii 

COS0 IQ 


. , , sin* 0 sin* <& , 

sin^ 0 COS'* + • - ■ ^ + cos^ 0 


(l-e^). 


> 


for all values of directivity parameters (N) and (e). It may be verified that: 

N = 22.8, € = 0.995827 


(26) 


i 

I 


(27) 


for an elliptical beam with -3 dB beamwidths of 19.88® and 9.0®. After Euler-angle rotation 
(a 90®), the semimajor axis of the elliptical beams is parallel to vector ly, as shown in Figure 14. 
The parameter N = 22.8 results in a -15.0 dB feed edge-taper in the y-sense for the central beam 
(i - 3) with 0 = 31®. The edge-taper in the x-sense was subsequently varied by controlling the exci- 
tation (S|) for the individual beams; this is identified for the individual radiation patterns presented 
in this document by means of a graphic symbol for both amplitude (Sj) and phase (i//.) of the 
beams (i = 1 , 2, . . ., 5). 
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Returning to Figure 13, it can be seen that the 9° beamwidth is associated with a corresponding 
feed dimension: 


70X 

D 

^ BW" (-3 dB) 


70X 

— 8X = 1.7 m 
9 


( 28 ) 


to the nearest integer multiple of half-wavelength. Defining far field for the feed in terms of 

Rc ** = 64X = 13 J m (29) 

it is concluded that the five prime-feed radiation patterns are everywhere Fraunhofer patterns, as 
seen from the surface of the offset paraboloid ( 7 ). It is noted that the radius of the Airy disk for 
- 8 ° plane-wave incidence, obtained by physical optics theory using a one-wavelength integration 
sampling interval and an annular disposition of samples over the refiector (y), is approximately 
1 .5 m. The inference is that the displaced Airy disk diameter (D^^ 3.0 m) is able to accommo- 

date the feed length (D'g. « 1 .7 m). 

Before the coma-correction simulation described above was implemented, several reference patterns 
were computed. The offset paraboloid wu first illuminated by a single pencil beam (i = 1 ), with 
-15 dB feed edge-taper, directed at the center of the refiector (C) (Figure 2). Feed directivity was 
controlled by parameters (N, e) = (22.84, 0.0) and beam-pointing by Euler angles (a, /3, 7 ) = (90.0°, 

3 1 .0°, 0.0°). The feed intensity was arbitrarily set at unity (Sj = 1 .0), and the initial phase at zero 
(^ j = 0.0°). A linear polarization moment (P(I) = 1 . 0 ) in the ( 1 ^^ ) direction was effectively con- 
verted to P(J) = 1 .0 in the (1^ ) direction under the assumed (a » 90.0°) Euler-angle rotation of the 
prime feed in this simulation. The feed displacement was taken to be zero = 0.0). 

A half-planar cut (0 = 1 80°) of the resulting far-field physical optics radiation pattern is shown in 
Figure 1 S. Beamwidth (degrees) is shown for -3.0, -10.0, and -20.0 dB relative to the beam maxi- 
mum. The first sidelobe was computed at -29.3 dB, and is due to the combination of feed taper, 
space divergence of the prime-feed fields, and the law of current induction on the siuface of the off- 
set paraboloid: 

K = 2 n X fi, . (30) 

This secondary beam is not squinted (B, =0°) and exhibits no aberrations since the feed phase cen- 
ter was taken to be at the focus (F) of the offset refiector under = 0.0. 

In order to verify the simulation when using five (i = S) elliptical prime-feed beams, a simple con- 
figuration was selected. The individual beam intensities were set to unity (S, = 1 .0) and the initial 
phases to zero (0, = 0.0°). Beam elliptidty was obtained as indicated by equation 27, (N, e) = 
(22.84, 0.995827). The central beam (i = 3) was directed at point C from the focus (F) by setting 
(fi^ - 0.0) and Euler angles (a^ ,^ 3 , 73 ) = (90.0°, 3 1 .0°, 0.0°). Effective polarization was again 
P(J) = 1.0. The incremental arigle (A/3) was set at a constant value (9.0°). 
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A half-planar cut (0 = 1 80“ ) of the resulting far-field physical optics radiation pattern is shown in 
Figure 1 6. Beamwidths are given at three power levels as before. Note that the half-power beam- 
width is slightly smaller than in Figure IS, which is to be expected since the individual beam inten- 
sities were equal. The first sidelobe level was -1 8.0 dB. This result agrees with the fact that uni- 
form illumination for a circular aperture is associated with a -1 7.6 dB first sidelobe level. The value 
e - 0.995827 has introduced a slight taper in the vicinity of the reflector edges (P = 9°, 53“). It was 
concluded that this mathematical simulation for five beams was correct and could be used to 
explore coma correction techniques. 

The prime feeds used to compute the radiation patterns of Figures IS and 16 were next translated 
to a point in the approximate focal plane of the offset paraboloid: = (1 .2720, 0.0, -0.7949) 

meters. The angle /3 was updated from 31.0“ to 34.0“ so that the single pencil beam under (1=1) 
and the central beam under (i = S) would again be directed to point (c) on the paraboloid. Figures 
17 and 18 show the far-field or Fraunhofer radiation patterns, and beamwidths are identified as 
before. The coma lobes for these simulations are at -20.0 and -1 2.8 dB, = 4.3“, and should be 
compared to -29.3 and -1 8.0 dB, = 0.0“. 

Prior to attempting coma correction, by introducing discrete phase differentials among the five 
elliptical beams, a conservative discrete amplitude taper was selected for these prime-feed beams: 

Sj = Sj = 0.316, S3 = 1 .0, Sj = = 0.562. This program input had the effect of producing a pro- 

gressive taper of -5 and -10 dB to either side of the central beam. The result is given by the radia- 
tion pattern of Figure 19, which shows a -19.0 dB coma lobe. 

The objective of reducing coma by compensating the cubic phase gradient with five primary beams 
having a common phase center (pe^ - constant, i = 1 , 2 , . . ., 5) was approached by attempting to 
optimize the phase of the beams corresponding to i = 3, 4, 5 while holding the amplitude taper 
S, = Sj , S3 = 1 , S3 = indicated above. In addition, phase variations that were arbitrary and 
departed from the theory were explored to further verify program behavior. A few of the results 
are presented here.. 

The phase parameter ) was varied over wide latitudes for the selected amplitude taper among the 
five beams, and it was concluded that a value = -13S.0“) was very close to optimum. An angular 
increment of = 9.0“ was retained, and the central beam angle of ^3 = 36“ was eventually used. 
Coma reduction was foimd to be rather insensitive to the set of angles (0.) since the latter affect 
illumination intensity, spillover, etc. for beams with a common phase center. A coma lobe of -22 
dB level was obtained by varying ) alone. (See Figure 20.) 

Similarly, the phases of (i//^ ) and (i/^^ ) were varied simultaneously, in accordance with and contrary 
to a cubic law across the aperture in a vertical (1 ^ ) sense. It was concluded that values = 
-130.0“, \jj^ - -30.0“ were close to optimum for coma-lobe reduction. A coma lobe of -28.3 dB 
level was obtained by varying and \l/^ simultaneously. (See Figure 21 .) It is noted that this 
beam exhibits approximately -30 dB first minima and sidelobe levels and beamwidth similar to 
those obtained in Figure IS, an unscanned case obtained with symmetric -IS.O dB feed taper. 
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The phases of , ^ 3 ) were varied siinultanebusly, and although some beam characteristics 

changed, no significant improvement over the result given as Figure 21 was obtained. At this stage 
of the exploration, a decision was made to calculate the orthogonal radiation pattern cut. Since the 
beam of Figure 21 is squinted (d, = 4.5°) an insitu orthogonal cut throu^ the beam maximum 
would imply both 9 and 0 variation with attendant plotting difficulties. Alternatively, a proper 
rotation of the entire discrete current distribution of the physical optics calciilation by 9^ = 4.5° 
would place the beam maximum on the (+1,) axis again. Consequently, any radiation pattern cut 
through the beam maximum can be made by setting 0 ^ constant and treating 9 as the single inde* 
pendent variable. 

The rotation of sheet current values (K) to obtain a new set (R'), 'as described above, was made 
simply by means of a qiiatemion transformation. A set of position vectors (fi) of (R) on the offset 
paraboloid was converted into a new set (p') bf (R') by the identical transformation. These opera- 
tions made possible the integration of sheet current, weighted by incremental reflector area, under 
the physical optics algorithm with minimal reprogramming effort. Vectors (R) and (p), above, may 
be regarded as “free” and “bound” vectors, respectively, in the present context. The same holds 
' true for their counterparts (K') wd (p'). 

A protracted discrission of quaternions is inappropriate here. It may be worth noting that a rota- 
tion in the amoimt 9^ in physical 3-space (E, ) appears as a half-angle or 9^ /2 operator in a trans- 
formation that resembles the similarity transformation of classical physics. The general rotation 
is executed about a preferred axis whose direction cosines are (oeq , /3q , 7q ), not to be confused 
with Euler angles (a, 0, y) used previously. In conclusion, the quaternion (Q^ ) simplifies to a vector 
(Wj ) of Ej , and 




•1 


(31) 


take 

Qj = + (iaQ+j0Q+k7Q)sin — (32a) 

and 

d 0 

Q\^ = cos - (i Oq + J0Q + k 7 q) sin . (32b) 

See References 17, 18, and 19 for background material. Note that (ij) equals (k), (k^ ) equals 
(-1), (kj) equals (-i), etc. 

The orthogonal pattern cut required here is made by setting = 9^, ~O,0q - 1-0, = 0.0, 

0 = 90°, 270°. Here, 9 is the independent variable as with unscaimed beams. 

Using the same program parameters as those for Figure 21 , the 8 = 90°, 270° cut effectively co- 
alesces the aperture distribution to a line parallel to the (1^ ) axis. A radiation pattern whose -10 
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and -20 dB beam widths are somewhat too broad was obtained by physical optics and is given as 
Figure 22. It shows a vestigial sidelobe of -22.0 dB level, about 6 dB worse than the corrected 
coma-lobe level and 3 dB better than the uncorrected coma-lobe level (Figures 19 and 20, respec- 
tively). The necessity of verifying final radiation patterns by means of principal plane and diag- 
onal plane cuts, at least, is apparent. 

ABERRATION CORRECTION (II) 

Attempts to minimize coma for increased plane-wave angle of arrival (J3 ^ -8°) at a frequency of 
1.414 GHz for the assumed geometry proved difficult. Since the effective aperture was diminished 
significantly, beamwidths increased. Estimation of for feed positioning by various means was 
tedious. (See Figures 10 and 11.) Previously, F' was taken to lie in the aperture plane (d 4°) 
and p, was determined accordingly. Hiis approach was now unsuccessful. A review of the assump- 
tions and an alternate approach were sought. The issue of apodization (taper) was set aside, and 
the phase aberration analysis was reviewed. 

The large wave-angle {d » -8**) indicated a radical departure from Gaussian optics, and the assump- 
tion of an offset paraboloid violated Hamilton’s condition that the instrument be symmetrical. It 
was noted that the aberration expansion of equation S provided for lateral feed displacement (h) 
relative to the system or focal axis (z), but exhibited no comparable axial feed displacement term 
(g). Also, the focal axis of the parent paraboloid and offset paraboloids did not coincide. Could 
(h') be assumed a lateral displacement with respect to z', the approximate focal axis of the offset 
paraboloid? Further, if the expansion S contained only particular orders of terms, was it capable 
of representing the phase aberrations in the aperture plane of the offset paraboloid for large wave 
angles (8 » -8®)? 

A ray-optics program was used to obtain the aberrations in the aperture plane of the antenna. The 
line cut (y = 0) in the ±x direction was selected from the aberration siuface, corresponding to 0 = 
0°, 1 80**, and a restricted expansion in powers of aperture radius (p) was made. It was tacitly 
assumed that terms up to and including the fourth power of p might suffice. Good closiu'e verified 


this assumption. In effect, the classical aberrations were rewritten: 

W,jjhpcos0 -^(W,„h)p (33) 

WojoP" (34) 

Wj 3 j hp’cos (/>■* (Wj 3 jh) (35) 

Wq^P-* -*Cffo4o)P* - (36) 


Preservation of the original concept of (h) orthogonal to z, or (h') orthogonal to z' hi the present 
case, was no longer required. The objective was simply to map coefficients, such as the associated 
terms (W^ 3 j h), in the focal region. This was done for both a parent and an offset paraboloid. 
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( 37 ) 


F/Dp = 18 ni/36 ra = 0.5 
F/D^ = 18m/15tn * 1.2 , (38) 

by selecting values of on a Cartesian grid in the (xz) plane. 

Figures 23 through 26 illustrate the coefficients associated with lateral focus, axial focus, coma, 
and spherical abenation of a symmetric parent paraboloid. It can be seen from Figure 23 that 
the classical aberration for lateral focus is observed along the system axis (z) for which (h) equals 
zero. For small nonzero values of (h), the lines 

Gj = 0^111 1») P ® constant (39) 

are nearly parallel to z. As the departure from Gaussian optics increases, (G, ) is seen to be con- 
stant along curved lines. Because the (cos 0) term has been suppressed, the coefficients assume 
positive or negative values. 

The axial focus coefficients are independent of lateral image displacement in the classical dis- 
cussions where axial image displacement (g) is not introduced. In the xz plane, 

Ip.l » (g*+h>)»4 . (40) 

The lines 


* constant (41) 

are shown in Figure 24. They assume positive or negative values also and are essentially circular 
arcs. 


The coma coefficients are displayed in Figure 2S. It can be seen that the classical aberration is 
observed along the system axis (z) for which (h) equals zero, as was the case for lateral focus. For 
small nonzero values of (h), the lines 

S4 * (Wj3j h)p® = constant (42) 

are nearly parallel to z, but become curves as the assiunption of Gaussian optics is increasingly 
violated. 

Finally, Figure 26 displays the spherical coefficients, also independent of (h) in the classical or 
receive sense. The lines for which 


Sj = (^040^ fi* - constant 


(43) 



appear to be parabolic arcs with positive or negative values. It may be verified that the contribu- 
tion of this aberration term, on an absolute basis, 

0 < p < 1 (44) 

. is rather negligible in the present discussion. 

Figures 27 through 30 illustrate the same coefficients as Figxues 23 through 26, but apply to the off- 
set paraboloid. Figure 27 suggests that the inclined focal axis (z') may replace z, and (h') may be 
taken orthogonal to z' in analogy with parent paraboloids, and this is also borne out by computed 
radiation pattern squint angles. Figure 28 shows essentially circular arcs and should be compared to 
Figure 24, but the center of those arcs is not along z'. Figure 29 shows curved lines somewhat 
similar to those in Figure 25, but symmetry about z' is lacking. Figure 30 departs drastically from 
the parabolic arcs of Figure 26, although some vestiges remain. Syimnetry about z' is lacking. 

In conclusion, it is difficult to obtain a close analogy between the parent and offset paraboloid be- 
havior. The asymmetric offset paraboloid has distinct featrues. Introduction of (h') and z' sets 
the lateral focus analysis on a one-to-one basis between the parent and offset paraboloids. The ad- 
vantage gained by mapping the modified aberration coefficients of Figures 27 through 30 is that a 
strategy for selecting (p^) can be developed and trade-offs between various aberration terms can be 
determined quantitatively. It is noted that the present approach was restricted to the line (±x, y = 0) 
in the aperture plane. The implication is that the aberration coefficients should, properly, be 
mapped for other lines (y 0) in the aperture plane. 

An'immediate simplification, predicted in Figure 27, is that the beam angle (d^) can be held con- 
stant while coma-lobe reduction is being carried out. It can also be seen that minimization of the 
term associated with coma, by intersecting the arc 

(W,3, h) = 0 (45) 


along a line 


(Wjjj h) = constant . (46) 

Figure 29 results in an enormous increase in the axial aberration. The departiue from the arc 

(Wo,o) = 0 (47) 

of Figure 28 is excessive. It may be instructive to locate the phase center of a multibeam array 
situated at 


Pj = (1 .68, 0.0, -2.0) meters 


(48) 
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in Figures 27 through 30 and compute the aberration values introduced. Once (p has been 
selected, the task of determining the relative phases (\]/^), between these beams with common phase 
center, may be undertaken to reduce coma. 

Figure 31 shows the result after optimizing the set of phase values for the feed position (p^) of 
equation 48. A coma lobe of -26 dB was obtained, 1 dB lower than the first local maximum at 
9 = 9.8° on the opposite side of the main beam. Beamwidths are given at four power levels. 

FOURIER TRANSFORMS 

As the present study of the correction of aberrations computed in the aperture plane of the offset 
paraboloid antenna progressed, the Fourier transform became ubiquitous. In anticipation of a more 
fully automated computer program in which the criteria of aberration correction are established in 
the aperture plane, rather than in the far*field region, the role of the Fourier trtmsform becomes 
increasingly important. Although two-dimensional transforms are required throughout, and even 
three-dimensional transforms could be used to advantage when dealing with the paraboloid, the sub- 
sequent discussion is restricted to one-dimension for convenience. 

Beginning with the relationship between the paraboloid antenna aperture distribution (AP^) and 
the far-field radiation pattern E (6, ^), the definition of the Fourier transform and its inverse are 
applicable as indicated previously (Ref. 16, p. 275 and Ref. 20, pp. 1 1 and 13): 


(49) 


(50) 


(51) 


and 


Alternatively, 


H(0 


<31 

I 




¥t) 


am 

£ 




if [h(t)l • H(0 


and 


[H(f)j « h(t) . (52) 

Next, as a direct consequence of these definitions, the Fourier transform isolates the Gaussian 
aberration, 


Wm h 


(53) 
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from other aberrations that affect image quality. Since the introduction of a linear phase gradient 
(e**"*®*) results in 


?(AP^) 


2F[h(t)e'""'“‘] = 



-j2 7T (f-fo)t 

f(t)edt =H(f-fo) 


(54) 


the shift theorem that ensures distortionless beam-swinging is established (Ref. 16, p. 282, and 
Ref. 21, p. 84). Equation 54 is sometimes referred to as the frequency-shifting property of the 
Fourier transform. 

Following this, the concept of a transform of a transform appears in several contexts. Given a 
prime-feed distribution (AP^) at focus (F) whose spatial Fourier transform leads to the paraboloidal 
antenna aperture distribution (AP^), to within a space divergence factor (r*^ ) and a minor perturba- 
tion due to the law of current induction (equation 30), it follows that: 

y[i(AP,)I = (AP^) = E id , «) . (55) 

Consecutive transforms may also occur with circuits and will be discussed in greater detail in con- 
nection with the Butler beam-forming matrix (Ref. 22 and Ref. 23, p. 258). Since the circuits may 
lead to the prime-feed distribution, cascaded transforms arise. A specialization of the concept is 
sometimes referred to as the symmetry property of the Fourier transform and follows from re- 
writing equation 50 as: 


h(-t) 



-j2irft 

H(f)e df 


and interchanging parameters t and f, so that 



-jairrt 

H(t)e dt = h(-0 


(56) 


(57) 


There is no tnmcation of H(t) in equation 57, which differs from equation 55. The concept of a 
transform of a transform is therefore also useful for extending Fourier transform tables without cal- 
culating new transforms (Ref. 21 , p. 85, Ref. 20, p. 32, and Ref. 24, p. 196). 


Of particular interest in the reduction of aperture aberrations are the Fourier transforms of uniform 
continuous and discrete distributions. The former is an idealization of the latter regarding prime 
feeds, but may also be used to generate a reference for paraboloids: 


H(0 = 



A e dt = 2AT^ 


sin (2n’T^f) 

(2jrT„0 


(58) 
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Here, aperture width equals 2T^, and A equals a real constant (Ref. 20, p. 24, and Ref. 21, p. 78). 
The above, partioxlarly relevant for the prime-feed array associated with the Butler matrix, is modi- 
fied by the shift theorem, and 





5(t-tje dt * e ‘ 




(59) 


Reference 21 , page 10, together with the linearity (superposition) property of Fourier transforms 
(Ref. 21, p. 82): 


■ gp [ajh, (t) + a,h, (t)] - a^Hj (0 + a^H, (0 . (60) 


Here 


gp [h, (t)] = Hj (f) 
and 


(61) 


SF [h, (t)l » H, (0 . (62) 

In the case of a uniform array of (n,) discrete sources a distance (t^) apart, with one source on the 
coordinate origin. 


Then 


and 


so that 




H(f) 


A(l + e * e 




H(0e ® = aU » + e 


■jaitf(at^) -jairf(n, 

**•►...•► e 




H(0 - » A 


(i 


(63) 


(64) 


(65) 


( 66 ) 
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From this result, 


H(0 = A 


_ sin[(27rft„)n,/2] 

1. sin[(2rfg/2] 


(67) 


prior to the application of the shift theorem, which introduces a linear phase gradient. This result 
for the linear array of discrete sources should be compared to equation 58, a continuous linear dis- 
tribution. If, in a physical application, each element width of the uniform array equals t^ , 

(2 T^) = (n, t^) = feed width . (68) 

Equation 6 may be generalized (Ref. 25, p. 451) for a cubic lattice array that is directly applicable 
in mathematical simulations of Butler array matrix feeds. 

Following Reference 25, let 


ka^ sin d COS 0 + 


(69) 

kaj Sind sin0 + 


(70) 

ka^ cos d ‘*’ 0 ^ 

9 

(71) 


where the a, are phase gradients and k is the wave number. The feed pattern for the cubic-lattice 
array is then: 


E>(0,«) = F^f,fjf3 


(72) 


where 


ann,7,/2 

anr,/2 


-l(ns-l V 




(73) 


and F^ is the source element pattern. Equation 73 may be compared to equation 67, which was 
obtained by Fourier transform theory. The cubic-lattice notation permits imaging in a ground plane 
(«3 = ir), but this capability is ordinarily defaulted (73 = 0 ) when the soiurce elements are electro- 
magnetic horns. A rectangular feed aperture of discrete sources becomes available by this analysis. 


Equation 72 may be regarded as a statement of the pattern multiplication principle or the result of 
time-domain convolution of a source distribution with an array factor distribution: 


lEy(9,0)l = (F„)(FjF 3 F 3 ) . (74) 
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The advantage of using equation 67 or its generalization (equation 73) instead of equation 25 is that 
there exists an inherent risk that the elliptical beams postixlated by certain parameters (Nj, e,, <3.) of 
equation 25 may not be realizable in the physical world. There is a reduced risk when using equa- 
tion 67, but uirrealistic specifications are still possible. For example, t^ may be selected so small that 
the phase centers of the discrete sources could never be brought adequately close together even with 
radical dielectric loading. For Butler array-matrix simulations, the element spacing is usually set at 
one-half wavelength. The size of the array and its distance from the surface being illuminated must 
still be regarded to justify use of the preceding far-f!eld radiation pattern equations for the prime 
feed in the simulation. 

BEAM-FORMING MATRICES 

The assumption of multiple beams having a common phase center has a basis in the physical world. 
Numerous beam-forming matrices exist (Ref. 23, Qiapter 3, p. 217). The Butler array matrix is one 
that forms (n^) orthogonal beams with a 4-dB crossover level, predicted on the array factor with 
half-wave separation between source elements. In systems simulation, it is convenient to express 
the properties of such a beam-forming device or circuit as a mathematical matrix. Scattering matrix 
notation may be adapted to transfer matrix notation, which admits chain matrix multiplication 
(Ref. 26, pp. 506-539). Computer processing for the circuit becomes a straightforward operation. 

A four-port 90° (square) hybrid junction is the basic building block of the Butler matrix, together 
with certain fixed phase shifters. Transfer matrix notation and the noncommutative order of multi- 
plication may be illustrated by a simple example for a lossless “weighting” circuit (Figure 32) whose 
output depends on a sin^e phase shifter (^j ). The example is given in more general terms, in- 
cluding a second phase shifter (^, ) and making provision for a transmission coefficient. This is 
done since cascaded Butler matrices have been proposed with variable phase shifters interposed 
(Ref. 27). Also, an attempt is made to explore the proposition that the Butler array matrix is a 
Fourier trarisform device. 

Figtire 32 shows the “weighting” circuit with inputs h, and h, , and outputs g^ and g 3 . The even- 
index quantities are absent because only the “transmit” submatrix of a partioned matrix is being 
considered here. Two null matrices and a “receive” matrix have been suppressed. The order of the 
matrices affecting the input signal is: . 


T,T3...T„.,T„ (75) 

so that the transfer matrix chain multiplication is carried out as: 

r = [T„] [T„.J ...{TJh . (76) 

For the “weighting” circuit, 
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-j//2 -1//2 

-1//2 -j//2 




o 


0 

„ • j '^2 
r,e 2 


H //2 - 1//2 
-l/Zz -j//2 




( 77 ) 


Defaulting the result, set equal to F, equal to unity, and equal to equal to zero. Then, 



3 

0 j 



s 




j 0 






(78) 


It can be seen that the magnitudes of outputs (g, and are the inputs h, and h^ , respectively). 
This is a symmetry aigiiment. Since a single square hybrid is an elementary Butler matrix, cascading 
two such hybrids might constitute “transform of a transform” in the Fourier sense. 


The proposition that a Butler matrix is a Fourier transform device is difficult to support formally 
for the two-port example shown here, or for the n^-port matrix, even though additional properties 
associated with a Fourier transform emerge. Suppose only a single hybrid is considered, so that: 


8t 

-1 

'j l' 


hi 

A 

“ /2 

J j_ 


^3. 


Suppose further that h, is identically zero, initially. Then it is possible to write 


and 




— 

/2 


*3 = 


-1 


/2 


/ 


(79) 


(80) 


(81) 


Now the single input (h^ ) equal to unity may be regarded as a delta- function and (ft^ ) may be taken 
to equal minus one-fourth for g^ and zero for g^ , so that: 

3P[«(t-g] = (82) 

reproduces the exponentials of equations 80 and 81 above. This has the appearance of aimexing the 
frequency-shifting property of Fourier transform theory to the symmetry property. It is noted, 
however, that the Butler matrix of n^ output ports accepts a single input, 5 (t - 1^^ ) and produces a 


20 



discrete set of output values e , 




-j2irf 


n 


s ° , whereas equation 82 indicates a 


continuous output. The discrete set of values may well be bounded by equation 82, but the prop- 
osition regarding the Butler array matrix as a Fourier transform device appears to be restricted. 


TRANSFER MATRICES 

Admitting the possibility of a restriction regarding the Butler array matrix as a Fourier transformer, 
the requisite signal processing may still be carried out conveniently using the preceding transfer 
matrix notation. In view of certain publications, including Reference 27, and the experience gained 
in the present effort, a transfer function consisting of three matrices would be useful in the develop- 
ment of future computer programs. The cascading of Butler matrices or power dividers of different 
dimension is anticipated. A single transfer function that accommodates these options may be 
written. Matrix conformability problems regarding Butler matrix dimensionality are circumvented 
by an artiGce. The higher dimension is carried throughout, certain parameters are defaulted, and 
the matrix of smaller dimension is thereby imbedded in the formulation. 

The circuit of Figure 33 cascades Butler matrices of dimension equal to eight and four. Notation is 
the same for the square hybrid example, previously given, as is the order of multiplication. An 
equivalent matrix for the entire transfer fxmction may be obtained by multiplication; 


[T„] » [Tj] [T,] [T,] (83) 

where [T^ ] is associated with a four^beam array matrix or power divider, [T, ] designates a bank of 
“weighted” phase shifters, and [T, ] is associated with an eight-beam array matrix. Prior to de- 
faulting any of the parameters. 
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(85) 


for the coma correction circuit. For the remaining inputs, assume that only hj is nonzero initially. 
It can be seen that product [Tj ] [li] results in a column vector (or matrix). The elements T 33 , 
T 43 , Tj 3 , and Tj 3 of [Tj ] are exponential terms of the stepwise linear phase gradient that is ordi- 
narily associated with an “antenna beam” from this Butler matrix. 
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Multiplying column vector [Tj ] [h] against matrix [Tj ] , in the order indicated, introduces ele- 
ments T33 , » Tj j , and of [Tj ] . The elements are weighted exponentials whose arguments 

are such that they remove the linear gradient due to [Tj ] and the comatic aberration as well. The 
weighting factors provide the amplitude taper of the elliptical beams of the previous discussion. 
Symmetric taper of the feed function alone is obtained with [T33 ] equal to [T^^ ] and [T^^ ] = 

[Tj j ] in matrix [Tj ] in the present discussion. Product [Tj ] [Tj ] [h] results in another column 
vector that, when multiplied against matrix [T3 ] , leads to four antenna beams. Each beam is asso- 
ciated with a column of matrix [T3 ] as indicated in equation 84 . 

CONCLUSION 

This document developed analytical/numeiical approaches to identifying and correcting the aberra- 
tions introduced by a general displacement of the feed from the focal point of a single offset para- 
boloid antenna used in radiometer systems. The transmit point of view led to an optical path 
difference in the aperture plane. Qassical aberration theory was then applied in an attempt to 
determine the mathematical composition of the optical path difference or total aberration. It was 
noted that the classical theory was useful, but large feed displacements and the inclusion of an 
asymmetrical antenna surface caused the investigation to depart from Gaussian optics and violate 
Hamilton’s assximption of a symmetrical instrument. No attempt was made to include the effects of 
apodizadon or amplitude taper in the geometric analysis here. 

Ray optics was also applied with the receive point of view, leading to the development of wavefronts 
and caustics, under the assumption of zero wavelength. This approach was augmented with the 
physical optics mapping of the focal-region fields at a wavelength of 0.212 meter to determine suit- 
able feed positions for a given plane-wave angle of arrival. 

Finally, physical optics was used for the transmit case to compute far-field radiation patterns and 
to test varioiis hypotheses. A common phase center was assumed for a prime feed that radiated five 
elliptical beams. Both phase and amplitude taper were used in a mathematical simulation that is 
also physically realizable. The prime feed was associated with the Butler array matrix and its con- 
straints. It was found that the quality of coma correction was set by feed size (number of discrete 
sectors on the reflector), and adverse effects orthogonal to the comatic phase gradient. Beam-scan 
an^e (projected aperture) determined minimum beamwidth. 

This document illustrates some of the methods that may be applied in future radiometric anteima 
design effort and points out limitations of those methods. The actual coma correction achieved in 
the various examples is considered incidental to the value of the methods discussed. Mapping of the 
aberration coefficients for parent and offset paraboloids has not been found in the open literature, 
and appears to provide useful insight for the coma-correction problem and other anteima studies. A 
modified aperture-plane analysis for complex distributions, combined with a Fast Fourier Trans- 
form (FFT) capability, could replace the physical-optics approach used in the present effort. The 
inclusion of transfer matrices is not imperative, but constitutes a useful alternative in an automated 
approach to coma correction. Finally, it appears that coma correction remains an art, rather than 
a science, since the computational path is very long. There are many variables and constraints; 
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came and effect cannot be expressed in any simple manner. An iterative approach appears to be 
necessary. 
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Figure 1. The optical path difference (OPD) for optical (receive )/microwave (transmit) views. 
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Figure 2. Geometry of the offset paraboloid. 
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Figure 7, Truncated analytical caustic curve (Tschirnhausen's cubic) for -4°, -8®, -16°, and -32 

plane-wave incidence (offset paraboloid). 
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Figure 11. Coordinates (p^ = (x^, y^, z^)l of displaced focus (F*) for -8.7* plane^ave incidence. 
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Figure 13. Offset paraboloid illumination of -8° plane-wave incidence 
showing feed (D,) and far-fteld radius (R^). 
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Figure 14. Illumination of offset paraboloid with five beams (i = 5). 
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Figure 15. Reference pattern. 






Figure 17. Reference pattern. 
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Figure 19. Secondary radiation pattern obtained with ampiitude-tapered five-beam primary feed. 
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Figure 21. Secondary radiation pattern obtained with phase/amplitude-tapered 
five-beam primary feed. 
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) for parent paraboloid. 
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APPENDIX A 


ABERRATIONS 


The optical invariants are: 


= h*, Ij » p*. I 3 = hpcos 0 (A-1) 

and the optical path difference is: 

OPD = 5;» ^ (A-2) 

i 

Let the Cartesian aperture coordinates be 

xapcos^, y»psin 0 (A-3) 

Then the classical “aberrations” of Gauss, Seidel, and Schwarzschild of mixed degree (2), (4), (6) 
are easily generated in terms of the optical invariants as follows. The aberration name is sometimes 
redundant in the literature.* 


*W. L. Wolfe, ‘The Infrared Handbook,** Environmental Research Institate of Michigan, Library of Congress Cat. Card No. 77* 

90786, 1978; and E. H. Unfoot, Recent Advances in Optia^ Oxford, 1955. 
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APPENDIX B 


PHYSICAL OPTICS SIMULATION (RECEPTION) 


E (x,y^) 




(nXH,)i//ds (B-1) 


H (xjr^) 


4n 



X 


V ^ ds 


iff S 


(B-2) 

(B-3) 



(B-4) 


Parameters: 

Frequency: f = 1.414 GHz 

Diameter of reflector: D„ = 15.0 m ”* <r,^„ = 7.5 m 
Central blockage: = 0.0 m 

Reflector displacement: (Xj^, y,p, 2 jp) “ (-10.5 m, 0.0 m, -18.0 m) 

Focal length: F« 18.0 m 

Source function amplitude, phase, polarization: (S, , ^ ^ , P^ (D) - (1 -0, 0.0°, 1 .0) 
Source directivity factors: (N, e) = (0.0, 0.0) 

Soiirce displacement: (pe) = (x^, y^, z^) = (-22.689 X 10® m, 0.0 m, 148.274 X 10® m) 
(Range = 1 50,000 meters) 

Source inclination (Euler angles): (a,/J, 7 ) = (90.0°, 0.0°, 351.3°) 

Integration sampling interval: (LI) - 1.0 (X), (Fig. 9) 

Source fimction: 




cos^ 9 

sin® 0 cos® + sin® 0 sin® 4 + cos® 0 ^ 


1.0 


under N - e = 0.0 for plane-wave simulation (remote spherical source). 
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